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CN" : ABSTRACT 

o 

■ The Wide-field Infrared Survey Explorer Preliminary Data Release Source 

^ ; Catalog contains over 257 million objects. We describe the method used to flag 

variable source candidates in the Catalog. Using a method based on the chi- 
^ ■ square of single-exposure flux measurements, we generated a variability flag for 

each object, and have identified almost 460,000 candidates sources that exhibit 
significant flux variability with greater than ~ la confidence. We discuss the 
flagging method in detail and describe its benefits and limitations. We also 
^ ■ present results from the flagging method, including example light curves of several 

types of variable sources including Algol-type echpsing binaries, RR Lyr, W UMa, 
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P ■ and a blazar candidate. 
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INTRODUCTION 



The Wide-field Infrared Survey Explorer (WISE; IWright et al. Il2010[ ). mapped the 
entire sky in four bands centered at wavelengths of 3.4, 4.6, 12, and 22 /im (hereafter Wl, 
- W2, W3, and W4). WISE conducted its survey using a 40 cm cryogenically cooled telescope 

O ' equipped with four 1024x1024 infrared array detectors that simultaneously imaged the same 

^ 47'x47' field-of-view on the sky. WISE flew in a 531 km sun- synchronous polar orbit and 

employed a freeze-frame scanning technique whereby the telescope scans the sky continually 
at a rate of approximately 3.8 arcmin/sec and a scan mirror freezes the sky on the focal 
Vh ■ plane detectors while 7.7 sec (Wl and W2) and 8.8 sec (W3 and W4) exposures are acquired. 

- - The FOV of each successive exposure set overlaps the previous one by 10%, and the scan 

paths on adjacent orbits overlap by approximately 90% on the ecliptic. The WISE survey 
strategy alternated stepping the scan path forward and backward on subsequent orbits in 
an asymmetric pattern that approximately matched the orbital precession rate. In this 
way, each point near the ecliptic was observed on every other orbit, approximately each 191 
minutes, and typically 12 independent exposures were accumulated for each point near the 



^Infrared Processing and Analysis Center, 
California Institute of Technology, 
Pasadena, CA 91125, USA; dlioffman@ipac.caltech.edu 



-2- 



ecliptic. The number of exposures increases with echptic latitude, reaching over 1000 at the 
echptic poles. 

The WISE observing cadence is well-suited for studying periodic variable stars with 
periods less than ~ 2 days near the ecliptic. The most common variables in this category 
are Algol, /3 Lyrae, and W Ursae Majoris (W UMa) -type eclipsing binaries, and RR Lyrae 
pulsating variable stars. W UMa variables are the most numerou s among thes e, as 95% 



of all variable stars in the solar neighborhood are W UMa stars (JEggen 1119671 ). but their 
mid-infrared emission is relatively weak, and their variation amplitudes can be small and 
difficult to detect. Therefore, Algol and RR Lyr variables are the most common periodic 
variables detected by WISE. 

Non-periodic variables are also detected by WISE. Among these are young stellar objects 
(YSOs) such as those seen in Orion and other star formation complexes, and active galactic 
nuclei (AGN). Blazars are be the most common type of A GN with variabilit y detected by 



WISE because of the short timescale variations they exhibit (jStein et al. Ill976l ). Other types 
of AGN vary in brightness over longer timescales, and therefore it is difficult to probe their 
variability with WISE, except perhaps near the ecliptic poles where the coverage period is 
longest. As we show below, AGN can be separated from the other types of variables using 
WISE colors in the cases where the source has a 4.6 /xm detection, and with more confidence 
when a 12 yum detection is available. 



The WISE Preliminary Data Release Source Catalog (ICutri et al. Il201ll ) contains en- 
tries for over 257 million objects that were detected on the coadded WISE Atlas Images 
covering 57% of the sky. Photometry for the Source Catalog entries was performed by 
fitting point spread functions (PSFs) simultaneously to all of the individual exposures cov- 
ering an object. The magnitudes reported for each object effectively average over any time 
variations, although measurements of the fiux variance across all individual exposures are 
available from which we can test for variability. Each source record is accompanied by a 
variability fiag, var_flg, that is related to the probability of fiux variation based on the 
distribution of source fiux measurements. 

In this paper, we describe the derivation of the WISE Source Catalog variability fiagging, 
and demonstrate how to apply the variability fiag to identify true variable sources. Examples 
of light curves for various classes of objects are presented, as well as some representative 
phased light curves for periodic variables with large numbers of measurements. 
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2. Variability Flagging Method 



The method used to identify the variables makes use of statistics computed in the multi- 
frame data processing. Traditional variability indices such as the Stetson L, J, and K indices 



(IWelch &: Stetson Ill993l : IStetson 1119871 . Il996l ) were not used. The main reason is that access 
to the single-exposure flux measurements is not practical on such a large scale, whereas 
statistics from the multi-frame generation step are readily available. The Stetson indices are 
also rather insensitive to transient events, such as flare stars, although the method we used 
also suffers from this problem. WISE was not designed to identify variables, so using the 
multi-frame statistics gives the most practical and useful variability information. 

During the WISE multi-frame processing, fluxes are measured from all individual ex- 
posures simultaneously at the location of the source detection in the coadded image. The 
standard deviation of the ensemble fluxes are computed and included in the WISE source 
catalog. From these data, a variability flag is generated {var_flg), which is an integer value 
(0-9) for each WISE band. A value of zero indicates insufficient or inadequate data to deter- 
mine a variability probability, while values of 1-9 indicate increasing likelihood of variability. 

As a first step, measurements are pre-filtered to eUminate sources that are not suited 
to make a variability determination. The following lists the pre-filter conditions used and 
the definitions and/or reasoning behind each. Each pre-filter was based on constraints on 
columns from the WISE Preliminary Release Source Catalog, with the exception of the 
single-frame SNR proxy. See Table [1] for a brief d escription of t he co lumns discussed in 



this paper. More detailed definitions can be found in lCutri et al. I ( 120111 ). In the conditions 



below and in Table [H "?" denotes a WISE band, and is equal to 1,2,3 or 4. 



> 0.8 : Reliability of single-frame measurements. This helps to ensure there is a 



detection of the source in the single-exposure image. These parameters ensure at least 
80% of the frames have a single-exposure detection at the location of the source in the 
multiframe image. However, this does not ensure there are also single-exposure fiux 
detections for the source, as the single-exposure and multi-frame fiuxes are measured 
at slightly different positions. 

w?rchi2 < 5.0 : This reduces confusion issues by rejecting sources with high 
template-fit chi-square values. High chi-square values normally indicate confusion with 
nearby objects, or that the source is extended, both of which are conditions that can 
lead to measurement contamination. 

cc_flags = in a least one band : Contamination by image artifacts, especially 
diffraction spike artifacts, can produce spurious variability since the artifact location 
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and intensity can vary between frames. Variability flagging is unreliable in the presence 
of artifacts and is avoided. 

• wlm > 10 : Small depth of coverage produces unreliable statistics since outliers have 
a larger effect on the overall variability. This filter helps to ensure that any variability is 
robust against noisy measurement outliers, by requiring a minimum depth of coverage. 

• Single-frame SNR proxy > 5 : This helps to ensure there is a single-frame detec- 
tion. See the more detailed explanation later in this section. 

• wlsnr > 5 : The coadd photometry SNR must be significant, otherwise there will 
be no single-exposure detections. This filter ignores sources that are too faint to have 
any single-exposure detections. 

• na = : Eliminates sources with active deblending, which indicates one more more 
close neighbors that can produce false variability. 

• nh < 3 : Limit to the number of blend components to two. This also avoids confusion, 
which can generate variability. Tests showed that some known variables were still 
detected with nh = 2, so the threshold was set at this value. For sources with nb > 2, 
spurious flux variations may be seen due to confusion. 

• volsigmpro not null in at least one band : There must be a multi-frame detection 
of the source. 

• wlsigpl not null in at least one band : There must be single-frame measurements 
for the source. 

The single-frame SNR proxy is needed to ensure the variability flag is applied to actual 
single-exposure source detections. The single-frame SNR cannot be estimated simply by 
using the value wlsnr/ywlm^ as the multi-frame SNR estimation employs a noise model 
based on zero-mean random fluctuations and a systematic component due to PSF error. The 
latter does not average out over multiple measurements, so scaling the noise variance by the 
inverse of the number of measurements is not applicable. To adjust for this, approximately 
100,000 random sources were used to compile the actual SNR of both single- and multi-frame 
measurements. A linear flt was then applied to the function 



Y(w?mpro) = vwlm 



wlsnTr 



where wlmpro is the coadd magnitude of the source for a given band, wlm is the number 
of frames that went into the multi-frame image, wlsriTse is the SNR of the single-exposure 
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measurement, and wlsnrc is the SNR of the coadded source. The result is that 

_ Y{wlmpro) ■ wisnrc , . 

ywlm 
gives a reasonable estimation of the single-frame SNR as a function of magnitude, given only 
the SNR of the coadd. 

An estimate of the population standard deviation, a, is determined in each band by 
taking the 65*^ percentile of wlsigpl in magnitude bins of 0.50. The 65**^ percentile was 
used instead of the median for conservatism and to ensure that candidate variable sources in 
the highest confidence level bins do not have many spurious variable sources mostly arising 
from unflagged artifacts. These distributions form a "null" non-variable distribution from 
which approximate which chi-squares can be calculated. These distributions are saved as a 
look-up table as functions of M and magnitude and later interpolated upon retrieval. There 
is an effective faint-magnitude cutoff at 16.0, 14.75, 10.75, and 7.0 for bands 1,2,3, and 4, 
respectively. This is due to the lack of usable data at these fainter magnitudes to make the 
look-up table, as these magnitudes are approximately the single-exposure detection limit. 
Sources fainter than these magnitudes will always receive a variability flag value of "0" in 
the appropriate band. At the bright end, there was also a general lack of usable data as 
sources start to saturate and have greater measurement errors as well as large chi-square 
values. The bright end of the look-up table ends at 6.5, 6.0, 2.5, and 1.5 magnitudes for 
bands 1, 2, 3, and 4, respectively. Linear interpolation was used for sources brighter than 
these limits, which imposes a higher than normal a for the calculation with these sources. 
This ensures that very bright sources are not flagged as variable in error due to sparse data 
and high measurement errors due to saturation. However, very bright sources flagged as 
variable should still be treated with suspicion. 

To generate the flag, we construct a function that is monotonically related to the likeli- 
hood that the null hypothesis is false, where the null hypothesis is that the observed fluxes 
came from a single real non- variable source. If the noise in the data were uncorrelated and 
Gaussian-distributed, the ideal function would be chi-square computed from the sample of 
fluxes. The chi-square statistic is given by 

X' = E ( ^ ) (3) 

where M is the number of measurements (equivalent to wim), Xi is the single-exposure 
flux measurement for frame i, /i is the flux population mean, and a^ is the flux variance 
of the non-variable population. This equation cannot be used as is, because the individual 
flux measurements are no longer available. The sample variance of the flux measurements, 
however, is available and given by wlsigpl: 
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M 

wlsigpl' = ^ (X, - ^,f (4) 



M 



Substituting wlsigpl'^ in Equation |3] and multiplying by w?m, the equation reduces to 

2 U7?m ■ [wlsigpiy 

The number of degrees of freedom, A^, is given by A^ = {wlm — 1). 

There are three significant approximations in Equation [5] that cause some deviation 
from the behavior of a true chi-square statistic. Appendix A expands on some of these 
approximations . 

First, it operates in the magnitude domain, whereas it is the flux domain that is best 
modeled as Gaussian. For sources with high SNR, the magnitude distribution is also well 
approximated as Gaussian, but for low SNR, the magnitude density function has an appre- 
ciable tail, positive skewness, and positive excess kurtosis (e.g., at SNR = 3, the skewness 
is 1.55, and the excess kurtosis is 3.63; at SNR = 10, these are down to 0.31 and 0.21, 
respectively). This causes the chi-square approximation to have excessive variance for low 
SNR sources. 

Second, in some cases the sample variance has been subjected to outlier trimming during 
the production processing. This appears to affect primarily bright sources and results in 
reconstructing the chi-square estimate with an overestimated number of degrees of freedom. 

Third, although correlated errors do not affect the expectation value of the sample vari- 
ances, ignoring them results in a significant inflation of the variance of the reconstructed 
chi-square defined in Equation 5, making it considerably less efficient than a true chi-square 
but still statistically a monotonic function of the likelihood of the null hypothesis and there- 
fore useful for our objective. Using it with the chi-square Q distribution propagates this 
approximation error, but as discussed below, the value of Q is presented with very low pre- 
cision, namely a one-digit value for — log^gQ, and so errors of a factor of 3 in Q operate at 
the truncation-error level. We require extreme statistical significance for variability tagging, 
so even if — logj^gQ is off by as much as —3, it is still useful as an indicator that the null 
hypothesis should be rejected. To calibrate an empirical Q function for the actual magnitude 
data is beyond the present scope, and so we use the chi-square Q function for this purpose. 

Because of the preponderance of faint sources, the dominant approximation error ap- 
pears to stem from the non-Gaussian magnitude distribution. Figure [1] shows the observed 
distribution of — logj^gQ as a function of chi-square for the case of 10 degrees of freedom and 
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Wl magnitudes of 13 ± 2. The observed distribution is very close to the simulated case for 
SNR = 5, which included the effects of the magnitude distribution being non-Gaussian but 
did not include the effects of ignoring correlations in the photometric errors. 

The probability, Q, that x^ would be at least as large as what was observed given the 
null hypothesis, is calculated from the integral of the probability density function: 



Q = l 



2^/2r(iV/2) 



^(iV/2-l)g-x/2^^ 



(6) 



where F is the gamma function. The value of the band component of the variability flag, F, 
is then defined as the floor of the negative logarithm of Q: 

F = floor [-logioQ] (7) 

with F clipped at 9 for F > 9 and set to 1 when F = 0. The final variability flag, var_flg, 
is a four character string comprised of the values of F for each of the four WISE bands. 

Figures I2E show the results of the method applied to a large subset of the WISE 
preliminary release data set. The figures contain wlsigpl as a function of magnitude for 
each of the four WISE bands. Each colored cross corresponds to a different value of F, with 
grey dots corresponding to F = 0. 

The values of F cannot be replicated using single-exposure flux measurements. The 
method uses the standard deviation of the flux measurements from single-exposures in the 
multi-frame pipeline. The fluxes are calculated at the position of the Catalog source, which 
can differ slightly from the position of the source in the single-exposure images. Therefore, 
the fluxes used in the calculation of wlsigpl are not the same as the fluxes in the single- 
exposure working database. 

There are some known limitations with the flagging algorithm. Confusion, close com- 
panions, and nebulosity can create false variability due to increased noise and positional 
errors compared to the single-exposure source locations. Some unflagged artifacts are still 
present in the preliminary release and contribute to contamination. The most common of 
these are sources falling within the halo radius or along a diffraction spike of a parent source 
that is off-image. The flagging method is not very sensitive to short-lived, low-amplitude 
transient events (e.g. flare stars and cataclysmic variables). A different method is required 
to detect these types of transients. Despite the largely successful attempt to estimate the 
single-frame SNR from the coadd SNR, magnitude, and number of frame overlaps, some 
sources are flagged as highly likely to be variable and have few, if any, single-exposure de- 
tections in the given band. The population distribution of wlsigpl"^ does not exactly follow 
a gaussian or chi-square distribution, but is close to a chi-square where the upper-percentile 



cut of 65% for 0" is a correction. And finally, artifacts are often flagged too aggressively, 
leading to many cases where F = where the light curve is actually quite usable. Many of 
these issues will be addressed in the final WISE data release. 



3. Results 

The WISE Preliminary Release Source Catalog consists of over 257 million objects, of 
which 459,906 are assigned var_flg values of 7 or greater (Eq. 5) in at least one band. 
Only about 5% of these objects have SIMBAD associations with known variable sources. 
The types of periodic variables recovered are dominated by RR Lyr, Algol-type echpsing 
binaries, and W UMa binaries. These types of variables are the most numerous due to the 
relatively short periods and high amplitudes. There are few variables seen in WISE bands 
3 and 4 mainly due to the lower sensitivity and the declining spectral energy distributions 
of many stellar sources at these wavelengths, but also because the reddest objects tend to 
have longer periods than can be observed with one WISE epoch (e.g. Miras and AGN). 

To extract the light curves of suspected variable sources, the position of each catalog 
entry is queried against the WISE single-exposure source database via the GATOR query 
service of the NASA/IPAC Infrared Science Archive (IRSA). Sources in the single-exposure 
database are not cross-matched with the multiframe sources, thus the association is done by 
a narrow position cone search. A search radius of 1.5 arcsec was used to generate the light 
curves in this paper. In crowded regions, or for sources with a close neighbor, it is possible 
for measurements of a non-target source to be included. In these cases, one must examine 
the cases where multiple sources are extracted from the same frame and decide which source 
is the target object. 

The typical periods of RR Lyr variables are between 0.2 and 1.0 day, thus more than 
one complete cycle is usually observed by WISE. An example light curve of an RR Lyr is 
given in Figure Ei There are many similar light curves available in the WISE data, but care 
must be taken to identify these properly as RR Lyr variables. A common problem is the 
similar light curve shapes of W UMa and RR Lyr variables of the RRc subclass. The RRc 
subclass objects are believed to be first overtone pulsators and have nearly sinusoidal light 
curves, which can be confused with other types variables. The RRab subclass (fundamental 
mode, asymmetric light curve) is more easily identified. RR Lyr variables are important 
distance indicat ors, and WISE can help ex tend the metallicity-period-luminosity relation 



into the mid-IR (jSandage fc Tamman II2006I ) 



Algol-type eclipsing binaries typically have longer periods ranging from fractions of a 
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day to many days, thus most WISE Algol light curves are partial, often consisting of just one 
eclipse. In these cases, period identification is impossible. However, for short period Algols 
and those near the ecliptic poles, where coverage duration is longer, period identification 
is possible. Figure [7] is an example of an object at high ecliptic latitude with coverage of 
a total of 4 eclipses. The object has not previously been identified as variable and is a 
candidate Algol-type eclipsing binary with WISE designation WISEP J092844.91- 732635.8 
and a period of 3 .1282 days, determined from chi-square phase dispersion minimization 



(]Stellingwerflll978[ ). 



The periods of W UMa binaries are typically between 0.25 and 1.2 days, so WISE often 
observed several complete cycles. Many observed W UMa binaries have enough cycles where 
a phased light curve is possible. Figure [H] is an example of such a case, where the W UMa 
binary V592 Car is phased to its orbital period of 0.8011 days. 

Sources with a high variability likelihood in both Wl and W3 as well as red colors in both 
W2— W3 and Wl— W2 are likely blazars or YSOs. Blazars and YSOs have distinguish ing red 



color and non-periodic variability in all bands, when available ( jWright et al. II2010I . Figure 
12). This is clearly seen in the color-color plot of Figure M The dots are sources with 
a variability flag of "9" in both Wl and W2, whereas the red symbols are sources with 
a variability flag of "9" in Wl and W3. All of the known blazars and YSOs are in this 
color regime, and this simple color cut when combined with the variability flag can quickly 
identify blazar and YSO candidates. The colors of blazars and YSOs can be quite similar, 
and distinguishing between the two can be difficult with WISE data alone. However, galactic 
latitude and observations at other wavelengths can generally be used to rule out a YSO. An 
example of such a candidate is shown in Figure [TOl The source has the WISE designation 
WISEP J053158. 61-482736.1 and has WISE colors and light curve structure consistent with 
a blazar or variable YSO. The position of t his source is very cl ose to that of the known 



flat-spectrum radio source PMN J053 1-4827 (JGriffith et al. Ill993l ). which suggests strongly 



that this object is blazar-like rather than a YSO. Variability in WISE bands 3 and 4 is rare, 
but tends to be dominated by these types of red, non-periodic variable objects. Correlation 
of candidate red WISE variables with radio source lists may provide an effective means of 
discriminating between the two prior to spectroscopic follow-up. 

It is difficult to determine the fraction of false variables, as the vast majority of the 
objects are unknown. However, approximately 20% of F = 7 sources appear to suffer from 
confusion noise and unflagged artifacts. Figure [11] is an example light curve of a source with 
var_flg=^7110\ which is likely a false-positive due to confusion with a nearby source. The 
false-positive rate decreases quickly with higher F values, with only about 2% of F = 9 
sources having the same problem. Most of the data in the preceding figures are included 



-lo- 



in the WISE Preliminary Release, with the phased light curve figures containing data not 
included in the preliminary release, but that will be in the WISE final data release. 



4. Summary 

We have described the method used to identify nearly 460,000 objects as significantly 
variable in at least one band in the WISE Preliminary Release Catalog. The method uses the 
ensemble statistics available during the multi-frame processing, and assigns a numerical value 
to each source and band that corresponds to the probability that the source is consistent 
with fiux variations seen in the bulk "non-variable" population within the same magnitude 
bin. While there are some limitations, the method is practical and effective, given the format 
and generation of the WISE data. For the WISE All-Sky Release, variability flagging has 
been further improved, incorporating band-to-band cross correlations, and contains over 6.8 
million objects that are signiflcantly variable. 

The objects identifled as signiflcantly variable are biased toward periodic variables, with 
RR Lyr and Algol-type variables being the most common. Light curves of variables with 
short periods and/or high depth of coverage can be phased, and period flnding or period 
reflnement is possible in these cases. WISE band 3 variables contain many AGN and YSO 
sources, which we show can be separated from many other types of variables by a simple 
color selection. WISE provides the best all-sky variability data in the mid-IR to date, and 
will serve the community well. 

This publication makes use of data products from the The Wide-field Infrared Sur- 
vey Explorer, which is a joint project of the University of California, Los Angeles, and 
the Jet Propulsion Laboratory/California Institute of Technology, funded by the National 
Aeronautics and Space Administration. Long-term archiving and access to the WISE single- 
exposure database is funded by NEOWISE, which is a project of the Jet Propulsion Labo- 
ratory/California Institute of Technology, funded by the Planetary Science Division of the 
National Aeronautics and Space Administration. This research has made use of the NASA/ 
IPAC Infrared Science Archive, which is operated by the Jet Propulsion Laboratory, Cal- 
ifornia Institute of Technology, under contract with the National Aeronautics and Space 
Administration. 

The authors would like to thank the anonymous referee for valuable comments leading 
to an improved evaluation of the approximation errors affecting the variability indicator. 
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A. Appendix: Sources of Approximation Error in the Variability Indicator 

Photometric errors are often weU described by Gaussian statistics, and this is the case 
for WISE fluxes measured in DN, i.e., data numbers proportional to flux. For this case, the 
probability density function describing the distribution of possible flux measurements of a 
source with true flux F is 



PfU) = ^^ (Al) 

Astronomical magnitude m is related to flux / according to 

m = -2.5 logio {j\ = ^0 - 2.5 log^ol/) (A2) 

where vtiq is a magnitude zero point corresponding to a reference flux /o that is in the 
same units as /. The second form of the relationship is the one most often encountered 
in practice, and so we will use that. Note that in taking the logarithm of the flux, we are 
implicitly clipping the domain to / > 0. 

Since m is a function of the random variable /, we can use standard theory of functions 
of random variables to obtain the probability density function for m. 



Tnn —m 
10 275 -F 



2 



, , ln(10)10(^^^)e ^f 
Vm{m) = ^ --j= (A3) 

For relatively high SNR, the shape of this density function is practically indistinguishable 
from a Gaussian, but for low SNR, it differs very signiflcantly, being positively skewed and 
leptokurtic. Since most sources detected in a survey are in the neighborhood of the detection 
threshold, they dominate the statistics, and their large deviations from Gaussian behavior 
cause the variability indicator based on a chi-square formula to have large tails that distort 
the Q distribution relative to a true chi-square. 

The variability indicator is based on a reconstruction of the expression 

v^ [mi — my , . ,. 

^^E \2 (A4) 



CJ 
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This expression is chi-square if and only if the rrii are Gaussian and uncorrelated. We 
denote this function \1/ instead of x^ because we have seen above that the rrii are not Gaussian 
for low SNR, and now we examine the error incurred by ignoring their correlations, which 
they develop to some extent because of the flux estimation algorithm. Here we consider only 
the latter effect and take the actual distribution to be Gaussian. That allows us to use the 
general multivariate density function for A^ correlated Gaussian random variables, the vector 

TUN 



-xir/2 



PArfrMjv) = : , (A5) 

where Q^ is the covariance matrix, and \^n\ is its determinant. Using pij to denote the 
correlation coefficient for rrii and rrij, the covariance matrix and Xn ^^^ 



Q 



N 



( 



cr; 



P2lO'2Cri 



Pl20'l0'2 

a' 



XPNlCTNCri PN2<^NCr2 



Pino-i(Jn\ 

P2N<^2<^N 



a 



N 



J 



(A6) 



X% 



N N 

y y Wij{mi — fn){mj — fn) 



where the Win are the elements of the inverse of the covariance matrix. 



(A7) 



For any value of A^, we can compute the mean and variance of \E' by formally performing 
the corresponding moment integrals using the joint density function above. The result is that 
the expectation value of \l/ is A^, the same as for a chi-square random variable, independently 
of any correlations, which is obvious from the fact that the definition of \E' consists of a sum 
of A^ terms, each of whose expectation is unity, and the variance of \E' is 



N-l N 



var(^) = 2N + 4 ^ ^ 



P§ 



(A8) 



i=i j=i+i 



If the correlations are all zero, then this variance reduces to 2A^, and ^ is chi-square. 
Otherwise, any correlation can only increase the variance, causing Q{^) to approach zero 
more slowly than the Q distribution for chi-square with correlations properly taken into 
account (e.g. Equation A7). 
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Fig. 1. — Cumulative distribution of the empirical — logj^Q Q values (stepped histogram) 
for a depth-of-coverage of M = 11 and the Wl magnitude range: W < m < 15, where 
Q = probability > x^- Also shown is the expected trend of — logj^gQ predicted from a x^ 
distribution for M — 1 = 10 degrees-of- freedom (dashed). The dotted and dash-dotted curves 
are from Monte-Carlo simulations to illustrate the impact of using magnitude units (instead 
of fluxes) on the noise-variance for two signal-to-noise ratios. The observed distribution was 
generated by a subset of the preliminary data release, consisting of approximately 6.6 million 
sources. See text for details. 
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Fig. 2. — The standard deviation of the single- frame measurements as a function of the Wl 
multi- frame magnitude (left panel). Each colored cross corresponds to a different F value. 
Grey dots are sources with F = 0, which lack sufficient data for a classification. The right 
panel is a histogram of the different F values with the appropriate bar color for all data in 
the left panel. 
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Fig. 3.— As with Figure El but for W2. 
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Fig. 4. — As with Figure |2l but for W3. 
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Fig. 5. — As with Figure El but for W4. 
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Fig. 6. — The light curve of the known RR Lyr-type pulsating variable RS Oct. 
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Fig. 7. — The phased hght curve of a WISE candidate Algol. The WISE designation is 
WISEP J092844.91- 732635.8 and is phased to a period of 3.1282 days, determined by chi- 
square phase dispersion minimization. 
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Fig. 8. — Phased WISE light curve of the known W UMa-type contact binary V592 Car. 
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Fig. 9. — Color-color plot of WISE selected variables. Filled black squares are sources with a 
variability flag of "9" in both bands 1 and 2. Red symbols are sources with a variability flag 
of "9" in band 1 and 3. YSO and blazar candidates can be identified by the red population 
far from the main (stellar) locus of sources. The black squares in this regime are also likely 
YSO/blazar candidates, but lack significant variability in W3. Known blazars are marked 
with a circled triangle, while known YSOs are marked with a filled red circle. 



-23- 



IN 






^ 
2 ^ 

55.0 



^ I E ■ 



I I I'^uiiii i|jiii|iniii "I 



.'i^v^',';^^ 



56.0 56.5 57.0 

MJD - 55200 



:5Z*5' 



58.0 



Fig. 10.— WISE light curve of WISEP J053158.61-482736.1 in all four bands. The non- 
periodic variable nature of the light curve, the galactic latitude, and the WISE colors are 
consistent with a blazar. The object has no previous classification. 
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Fig. 11. — Light curve of a source that is not hkely variable. The source has a variabihty 
flag of '7110', and has several close neighbors near its position, which likely affect the PSF 
photometry. 
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Table 1. Selected WISE Catalog Column Definitions 



Column Name 



Definition 



w?mpro 



The profile-fit photometetric measurement in magnitudes for the 
specified band. 



w?sigmpro 



The profile-fit photometric measurement uncertainty in mag units 
for the specified band. 



w?snr 



w?rchi2 



The profile-fit measurement signal-to-noise ratio. 

Reduced x^ of the specified band's profile-fit photometry measure- 
ment. 



w?sigpl 



Standard deviation of the population of the specified band's fluxes 
measured on the individual frames covering the source, in magni- 
tudes. 



w?m 



The number of usable frames containing the source's coordinates in 
the specified band. 



w?nm 



ccJlags 



The number of detections with SNR > 3 in the specified band. 

A four character string, one character per band [W1/W2/W3/W4], 
that indicates that the photometry and/or position measurements 
of a source may be contaminated or biased due to proximity to an 
image artifact. Non-zero characters denote possible artifacts. 



na 



Active deblending flag. "0" indicates the source is not actively 
deblended while "1" indicated the source is actively deblended. 



nb 



Number of PSF components used simultaneously in the profile- 
fitting for this source. 



varJlg 



The variability flag is a four-character string, one character per 
band, in which the value for each band is related to the probability 
that the source flux measured on the individual WISE exposures 
was not constant with time. A value of zero indicates insufficient or 
inadequate data to determine a variability probability, while values 
of 1-9 indicate increasing likelihood of variability. 



